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Q\ | Abstract 

OS ■ 

*"*> ■ Magnetohydrodynamic (MHD) waves are analysed in the early 

q-i Universe, in the inflationary era, assuming the Universe to be filled 

»JU ' with a nonviscous fluid of the Zel'dovich type (p = p) in a metric of 

&JQ' the de Sitter form. A spatially uniform, time dependent, magnetic 

field Bo is assumed to be present. The Einstein equations are first 

solved to give the time dependence of the scale factor, assuming that 



X 



j^ ■ the matter density, but not the magnetic field, contribute as source 



C3 



terms. The various modes are thereafter analysed; they turn out to be 
essentially of the same kind as those encountered in conventional non- 
gravitational MHD, although the longitudinal magnetosonic wave is 
not interpretable as a physical energy-transporting wave as the group 
velocity becomes superluminal. We determine the phase speed of the 
various modes; they turn out to be scale factor independent. The 
Alfven velocity of the transverse magnetohydrodynamic wave becomes 
extremely small in the inflationary era, showing that the wave is in 
practice 'frozen in'. 
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1 Introduction 

There exists the possibility that a primordial magnetic field was created at 
some early stage in the early Universe. This topic has attracted considerable 
attention in the recent past. The important point is that such a field may 
have left an observable imprint in galaxies. Our Galaxy, as well as many 
other spiral galaxies, is endowed with coherent magnetic fields, ordered on 
scales larger than about 10 kpc, with a typical strength of 3 fiG (= 3 x 10 -6 
gauss) |rj, [Q. This corresponds to a magnetic field energy density B 2 /8n 
which is of the same order of magnitude as the observed energy density p 7 
of the microwave background: p 7 ~ (4:fiG) 2 /8n [||]. 

There are two plausible explanations for the existence of these large scale 
galactic fields. The most popular explanation is that such fields are the 
result of a dynamo amplification || of a weak seed field ||, created some 
time in the early Universe, for instance at the electroweak phase transition 
0' 01- This transition took place at the instant t ~ 1CT 10 s, corresponding 
to a temperature of T c ~ 10 15 K (100 GeV). The size of the event horizon 
was at this instant H~J, ~ 10 cm. If the galactic dynamo is efficient at 
amplification, one estimates the magnetic seed field B seed to be lying in the 
range between 10~ 23 G and 10~ 19 G0, |Q. Another plausible explanation for 
the existence of the present magnetic fields is that they could have originated 
from a relatively large primordial seed field amplified by the collapse of a 
galaxy [|. 

Generally speaking, the magnetohydrodynamic theory (MHD) in curved 
spacetime is a relatively new development in astrophysics. We may mention, 
therefore, that useful recent review articles are given by Olinto || and Enqvist 
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Whereas in previous works it has often been assumed that the seed mag- 
netic field is created in the electroweak transition region, we will in the present 
paper go further back in time and consider instead the inflationary era. As 
is commonly assumed, this era took place from t ~ 10~ 35 s to t ~ 10 -33 s. 
The Universe was then in a state of violent expansion, with a cosmological 
scale factor being determined essentially from the cosmological constant A. 
If the influence from matter was negligible, as one often assumes, the met- 
ric was a pure de Sitter metric with scale factor proportional to exp ( J ^ t). 
During the inflationary period the Universe was subject to a large increase 
in size - one often assumes an increase of the order 10 50 - but the expansion 



might actually have been much larger. The temperature correspondingly fell 
from T ~ 10 27 K to T ~ 10 22 K. The Universe consisted during this era of 
a vacuum "fluid" with energy density p vac = A/8ttG, and extreme tensile 
stress, p vac = —A/8nG (this negative pressure giving rise to the repulsive 
gravitation). 

Brandenburg et al. [|l(J investigated, on the basis of a flat FRW fluid 
model, whether a primordial magnetic field decreases or increases with time. 



This work was related to prior work of Gailis et al. |JL1|1 . In the following, 
we will instead assume that there was a fluid of the Zel'dovich type in the 
early Universe. This means that the equation of state of the fluid was p = p, 
corresponding to a velocity of sound being equal to the velocity of light. 
There are reasons to expect that this extreme kind of fluid was present in 



the beginning, and at the end, of the inflationary era [0. Here, both p 



and p are taken to act as sources in the Einstein equations. In addition, 
we assume that there was a primordial magnetic field B present. Since the 
magnitude of Bo was so small, the magnetic field energy was too small to 
have any appreciable influence on the metric. We therefore examine simply 
MHD effects on a fixed metric, the metric being determined by A, p, and p. 
Note that because of the matter, the scale factor does no longer have the 
simple form exp (J^t) mentioned above. 

Our main purpose in the following is to establish the governing equations 
of the combined system, matter plus field. As far as we know, this has not 
been done before. We will establish the plane, linear, MHD wave modes on 
the spatially uniform de Sitter background. A characteristic property of this 
kind of system is that pR 6 = const. (R being the scale factor), instead of 
the conventional relationship pi? 4 =const. found in the radiation dominated 
FRW case. If B is perpendicular to k, we find that there is a longitudinal 
magnetosonic wave whose phase velocity is superluminal. Since this wave is 
nondispersive, its group velocity becomes accordingly also superluminal, so 
that this wave is not physically acceptable as an energy-transporting wave. 
Next, if Bo is parallel to k, there is a simple luminal pressure wave, and 
there is also a transverse magnetohydrodynamic wave propagating with the 
characteristic de Sitter variant of the Alfven velocity given by Eq. (50) below. 
In the inflationary era, the magnitude of the Alfven velocity is extremely 
small. 

Numerical investigations of the time development of a primordial mag- 



netic field, along the lines of flQ |, will not be undertaken in this paper. 

We adopt in the following Heaviside - Lorentz, instead of Gaussian, elec- 
tromagnetic units. 

2 The de Sitter - Zel'dovich scale factor 

We use the convention in which the Minkowski metric is r]^ = ( — h ++), 
Greek indices are summed from to 3, Latin indices are summed from 1 to 
3. We let f/ M = (U°, U l ) designate the four-velocity of the cosmic fluid. 

When the cosmological constant A(> 0) is included, it is convenient to 
write Einstein's equations on the form 

Gfiv = Rfiu — -^Rg^u = SnGT^, (1) 

where T^ is the modified energy-momentum tensor 

with T^ v = (p + p)U^U u +pg^u being the ordinary energy- momentum tensor 
of an ideal (non-viscous) fluid. Defining the modified energy density p and 
pressure p by 

~ P = P+ ^ P = p -^G> (3) 

we can write (2) as 

T flv = (p + p)U li U v +pg lu ,. (4) 

Consider now the line element having the de Sitter form: 

ds 2 = -dt 2 + R 2 (t)[dr 2 + r 2 (d6 2 + sin 2 6d<p 2 )]. (5) 

This form refers to comoving coordinates ( r is constant for a fixed mat- 
ter element). We consider the Einstein equations in the orthonormal frame 
(designated by "hats"). We need the following components of the Einstein 
tensor: 

G ^-§' *> = -T-W (6) 



The "energy" equation G° = 87rGT® now becomes 

U (7) 



o 
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whereas the "pressure" equation G\ = 8nGT~ becomes 



r r 
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- + w = ^Gp. (8) 

These are the field equations. We have also the energy conservation equation 
which for the de Sitter - Zel'dovich fluid leads to the behaviour p oc R~ 6 , 
instead of the behaviour p oc R~ 4 that is characteristic for the radiation 
dominated era. 

Let us now determine R(t), by first making the substitution 

p ( R-- & 



P0 V Kq J 

Here Rq and po refer to an initial instant t = to, which marks the beginning 
of the inflationary period. The first order equation (7) becomes 

S^(lf+f(l) 2 

Solving this equation for dt and integrating from t — to to t, corresponding 
to an increase of the scale factor from R to R, we get 



Ro 



/3A(t-t ) , S7rGp /A 



e vsA { t-t ) + ru/ sinh[V3A(t - t c 

1 + 8nGp /A 



:n) 



This expression gives the de Sitter - Zel'dovich scale factor. At first sight 
it may appear as if the Zel'dovich state equation p = p plays no role in 
Eq. (11) at all, since in the derivation of this equation we integrated the 
"energy" equation (7) which in itself makes no reference to the fluid's state 
equation. However, in the derivation we made explicit use of Eq. (9), which 
is characteristic for a Zel'dovich fluid. Therefore, Eq. (11) implicitly relies 
upon the Zel'dovich property. 

We see that for empty space (p = 0) the expression reduces to R(t)/R = 

exp[w-(t — to)], which is the conventional de Sitter form. For t = to, R(to) = 
Ro, as it should. 



3 Magnet ohydro dynamics in the early uni- 
verse 



3.1 Maxwell's equations 

As an introductory step, before embarking upon the energy-momentum con- 
servation problem it is instructive to consider Maxwell's equations them- 
selves, formulated in a curvilinear basis which we will choose in the present 
subsection to be the spherical basis. 

Let the spatial coordinates be numerated as x 1 = (r, 9, (p), corresponding 
to the basis vectors (e r , eg, e v ). The determinant of the spatial metric is, 
according to (5), 

7 = det( gij ) = R 6 r 4 sin 2 #. (12) 

We introduce the electromagnetic field tensor |T3|], [|TJ], JT5J, |T5| 



[in 



( -E r 
E r 

Eg ~^fB^ 



and the dual tensor density T> 



fW 



-Eg 



-Vl Br 



—Ey 
-V^B 6 

V^ Br 
o 



\ 



(13) 
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—By B r 

Bg -B r J 



(14) 



Three-dimensional indices are raised and lowered by means of the antisym- 
metric pseudotensor of rank 3: 



eijfc = 7 



V 2 6. 



ijki 



Jjk 



7~ l7 %fc, (15) 



m 



1. As ^oo = -1, 



-1/2 



5ijk being the antisymmetric Levi-Civita symbol with 5 
the gravitational permittivity and permeability are e = fi = {—goo) 
and we do not have to distinguish between the electric field E and the electric 
induction D, nor between the magnetic field H and the magnetic induction 
B. 



With J^ = (p, J) being the electromagnetic four-current density, the four- 
dimensional Maxwell equations 

iWi = o, (-g)- 1/2 ^\, = J" (16) 

can be expressed on three-dimensional form as 

(V x E) i = -^-d Q {^B% V-B = 0, (17) 

(VxB^ni^n V-E = p e , (18) 

the curl and divergence operators being defined as 

(V x E)' = e^djEk, V • B = 7 ~ 1/2 ^(V7 'B*). (19) 



3.2 Energy- momentum balance equations 

After having given above the general scheme for Maxwell's equations in curvi- 
linear coordinates, we shall in the following simplify the formalism by using 
"Cartesian" coordinates x,y,z, whereby the position four- vector becomes 
expressible as x M = (t, x, y, z). The line element is written as 

ds 2 = -dt 2 + R 2 (t)(dx 2 + dy 2 + dz 2 ). (20) 

The Christoffel symbols become simple: r°- = RR5ij, rL = r* = (R/R)Sij. 

Since goo = — 1 and g^ = the Lorentz factor takes the same form as in 
special relativity (cf. Sect. 10.3 in p~3fl): 
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;i-« 2 )- 1/2 , (2i) 



where u % = dx 1 /dt are the curvilinear components of the three-dimensional 
coordinate velocity u whose square is u 2 = gijU l uK The four- velocity of 
the fluid is U^ = 7(1, u l ). The energy-momentum balance equations for the 
mechanical subsystem are T" = f^, where T^ v is the energy- momentum 
tensor of the fluid as given above, and / M is the electromagnetic four-force 
density: 

U = V = (" E ■ J ' PeEi + (J x B),,), (22) 



the vector product in curvilinear coordinates being defined as 

(J x B)i = e ijk J j B k . 



(23) 



The energy-momentum balance equations for the fluid mechanical subsystem 
are 

U- (24) 



9 z 



We assume the fluid to be an ideal MHD fluid, implying that the electrical 
conductivity a — ► oo. It corresponds to a zero electric field E' = in a 
frame of reference moving with the fluid volume element considered. As J 
is finite, we then have that E = — u x B, and the volume force density f 
becomes equal to J x B. For \x = we get from (24), with the definition 
S = (p + p)7 2 u, 



—d^ip + p^-dop 



V-S + ^S-u 
R 



E- J. 



For fj, = i we obtain analogously from (24) 
1 



i? 3 



d (R 3 S i ) + d k (S t u k ) + d lV = (J x B)i. 



(25) 



(26) 



Consider next the electromagnetic subsystem. In the present coordinate 
system the field tensor schemes analogous to (13) and (14) are 



ji,U 



( ° 


-E l 


—E 2 


-E 3 \ 


E l 





R 3 B 3 


-R 3 B 2 


E 2 


-R^B 3 





R 3 B 1 


\E 3 


R 3 B 2 


-R 3 B 1 


/ 



(27) 



T- 



jW 



( R 3 E l R 3 E 2 R 3 E 3 \ 

-R 3 E l B 3 -B 2 

-R 3 E 2 -B 3 B l 

V -R 3 E 3 B 2 -B x / 



(28) 



With S^ v denoting the electromagnetic energy-momentum tensor the bal- 
ance equations for this subsystem can be written as —S v = / M , where 



s : 



F F va o v E aF a!i 



(29) 



In the present case where only the magnetic field contributes we have 

fi = (B k d k )B t - l -d t B 2 . (30) 

So far, the equation of state for the fluid has not been used. We now insert 
the Zel'dovich equation p = p, and assume small bulk velocities whereby 
7 2 PS 1. The term containing S • u in the \x — equation (25) is of order u 2 
and is negligible. We obtain the approximate equation 

^ (i? 6 p) + 2V-(p U ) = E-J. (31) 

The p = i equation (26) becomes analogously 

^d (R 3 pu t ) + d t p = {B k d k )B t - U t B 2 . (32) 

When expressed on vector form, 

-^d (R 3 pu) + Vp = (B • V)B - l -VB 2 , (33) 

this equation holds in an orthonormal basis as well. 

Equations (31) and (32) are starting points for the perturbative theory 
below. In addition we shall need also the first of Maxwell's equations (17), 
written as a vector equation. Observing again that E = — u x B, we have 

V x (u x B) = l-d (R 3 B). (34) 

R a 

It should be noted that when writing Maxwell's equations, we have not 
restricted the magnetic field to be weak. As far as the magnetic field is 
concerned, the only restriction made on it is that it is not so strong that it 
contributes as a source term in Einstein's equations. Our remaining physical 
assumptions, to summarize, are the presence of an infinite electrical conduc- 
tivity, nonrelativistic fluid velocities, and the Zel'dovich state equation. 



4 Magnet ohydro dynamic waves 

We will now consider plane, linear MHD waves on the spatially uniform 
background in the inflationary universe. The background is specified by the 
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magnetic field B , fluid density po, and pressure po- These quantities are 
time dependent. We shall need the proportionalities 

PoW«^,, B5(«)«^. (35) 

The first of these proportionalities follows from Eq. (9). The second follows, 
in the zeroth order approximation, from Eq. (34) when both sides of this 
equation are expanded around the background. Generally, we expand as 
follows: 

p = po + $P, B t = B l + 5B\ (36) 

with Sp and SB 1 being small quantities. 

From Maxwell's equation V x B = J + R~ 3 d (R 3 ~E) it follows that 

J = Vxffi + ^o[ J R 3 (uxB )] ) (37) 

which shows that J is a first order quantity. As E is also of first order, it 
is seen that the term E • J on the right hand side of (31) is negligible. We 
obtain from (31), (32), and (34) 

■^d (R 6 5p)+2p V-u = 0, (38) 

6R 

Zpout + -^PoUi + d t 5p = (B • V)<5B, - B • d^B, (39) 

V x (u x B ) = ^d (R 3 SB) (40) 

(recall that we are using curvilinear coordinates for which the curl and vector 
product operators are as given by (19) and (23)). 
We now expand the wave quantities as follows: 

M * = Mo e* e , Sp=^, 5B* = ^°, (41) 

where G is the phase 

Q = k^ = kiX 1 + k t = k-r-ujt, (42) 
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and u l , C, b l are constants. (We thus let the frequency u refer to the coor- 
dinate time t, not to the conformal time.) Substitution into (38) and (39) 
yields 

2p i? 6 k • u = Cu, (43) 

CR \ C 

1 + 3t R^ ) U °* ~~ B? h = (b ° ' Bo)/Cl " (k ' Bo)6 ° 1 ' (44) 

Since we are here considering waves moving on a "slowly" varying back- 
ground (in complete analogy with surface waves travelling on slowly varying 
currents or slowly varying topology in ordinary hydrodynamics), we can as- 
sume that w is much larger than the time derivative of the logarithmic scale 
factor: uo ^> (d/dt) InR. The second term in the first bracket in (44) is thus 
negligible, and we get 

C 

2ujp R 3 u 0i - —hi = (b • B )ki - (k • B )b 0i . (45) 

The third governing equation follows from Maxwell's equation (40): 

(k • BoK = (k • uo)SJ - |$. (46) 

We now take the vector k to lie along the x axis, and distinguish between 
two cases: 

Let Bo be perpendicular to k. From Eq. (45) it follows that W02 = «03 — 0, 
so that the fluid velocity uo is a purely longitudinal vector. From (46) it 
then follows that bo is transverse, collinear with Bo. The condition that the 
system determinant of the governing equations be equal to zero, yields the 
dispersion relation 

^ t 1 + S k2 (47) 

(k 2 = k^ki). From expression (42) for the phase it follows that the con- 
travariant component of the longitudinal phase velocity is v] = dx 1 /dt = 
uo jk\. We obtain 



v 1 l 

u long 



l + 5»!. (48) 



R\ 2p 
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The longitudinal phase speed is thus 



Vlong - >/9uV l(m g 



B„ 2 



N 1 + ^' (49) 



This expression is independent of the scale factor R. It has the same struc- 
ture as the standard expression vi ong = Js 2 + v\ for the longitudinal mag- 



netosonic wave in nongravitational MHD [0, s denoting the sound velocity, 
if we define the Alfven velocity as 

v A = Bo/ yj2^. (50) 

In the present case, s = 1. 

The following property of the dispersion relation (47) is however striking: 
the expression is nondispersive, which implies that the group velocity duo/dk 
in this case becomes superluminal. This is physically non-acceptable. We 
thus have to conclude that the Zel'dovich fluid model, being a maximally rigid 
model, does not permit the longitudinal magnetosonic wave to propagate as 
a physical wave, transporting energy. This result is not so unexpected, after 
all. The Zel'dovich model is very simple, and has in practice to be modified 
so as to take into account dispersive effects. 

Next, let Bo be parallel to k. There are two types of wave motion possible 
in this case. There is an ordinary longitudinal wave («02 = M 03 = 0, bo = 0) 
with contravariant phase velocity component equal to 1/R and phase speed 
accordingly equal to 1. There is also a transverse wave (uqi = 0), where b 
is antiparallel (if k • B > 0) to the transverse velocity u : 

b = - — (k-Bo)uo. (51) 

The dispersion relation becomes in this case 

^ = ^(k-Bo) 2 , (52) 

so that the contravariant component of the transverse phase velocity becomes 

vl = -t=Bl (53) 
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The speed of the transverse wave is thus 



Vtr = ~m=) ( 54 ) 



again an expression that is independent of the scale factor. We see that Eq. 
(54) actually agrees with the expression (50) for the Alfven velocity; as noted 
above, this correspondence is necessary in order to preserve the same formal 



relationship Vi ong = J s 2 + v\ as in nonrelativistic theory. 

To obtain a feeling of the order of magnitudes involved, let us estimate 
the value of po to be inserted in the expression (50) for the Alfven velocity. 
Recall that p is the density of matter and radiation; it does not include the 
vacuum energy p vac = A/SnG. We estimate po by considering the instant 
t = t\ = 1.4 x 10~ 33 s just after the termination of the inflationary era, when 
the Universe leaves the de Sitter phase and returns to the standard FRW 
radiation dominated form. When k — 0, the density is equal to the critical 
density 

Per = q T7- (55) 

(It is noteworthy that this critical density follows from Einstein's equations 
directly, for any metric that is expressible in the form (5); the expression does 
not rely upon a specific state equation.) Thus p = p cr sAt — ti. Taking the 
Hubble factor H = R/R during the inflationary era to be H = 5 x 10 34 s _1 
[Q, we obtain from Eq. (55) p cr = 4.5 x 10 75 g/cm 3 . This is an enormous 
mass density. The magnitude of the Alfven velocity (50) becomes accordingly 
extremely small: 

v A = 3x 10- 37 B (tesla) m/s = 3 x lO" 31 ^ (gauss) cm/s. (56) 

The transverse Alfven wave (a purely magnetohydrodynamic phenomenon) 
thus does not propagate appreciably under these extraordinary circumstances, 
unless the value of Bo is extremely high. And, as far as we know, this is not 
the case in the early Universe. 

5 Summary and final remarks 

We may summarize as follows: 
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(1) Our basic assumption is the presence of a Zel'dovich fluid, satisfying 
p = p, in a de Sitter universe whose metric is given by Eq. (5). The scale 
factor is given by Eq. (11). There is a uniform background magnetic field 
B present. The magnetic field energy is taken to be not so high that it con- 
tributes appreciably as a source term in Einstein's equations. The electrical 
conductivity is assumed to be infinite. The fluid velocities in the comoving 
reference system are taken to be nonrelativistic. The most actual example 
of application is the inflationary era, ranging from to — 10~ 35 s to t\ — 1CT 33 
s after big bang. 

(2) For a de Sitter universe, energy conservation considerations lead to 
the property pR 6 = const., instead of the conventional property pi? 4 =const. 
holding for a radiation dominated FRW universe. Considering plane, linear, 
MHD waves on the uniform de Sitter background, we find that the longitudi- 
nal magnetosonic wave (B _L k) implies a superluminal group velocity and 
accordingly does not correspond to a physical, energy-transporting, wave. Of 
the two remaining fundamental modes (Bo || k), the ordinary longitudinal 
wave (with vanishing perturbed magnetic field) can exist as a luminal wave, 
and there is a third kind of wave of magnetohydrodynamic type propagating 
with the characteristic Alfven velocity given by Eq. (50). 

(3) Estimating the value of po to be equal to the critical mass density 
of a flat FRW universe, Eq. (55), we find that the Alfven velocity becomes 
extremely small, Eq. (56). The magnetohydrodynamic wave in practice does 
not propagate; it is 'frozen' in the fluid. 

(4) The following remark ought to be made, concerning the physical mean- 
ing of the present theory. Our adopted state equation, p = p, refers to the 
content of matter, and radiation, in the early Universe. One may ask: is not 
the early Universe after all dominated by vacuum effects, corresponding to 
the state equation p vac = — p mc ? The point here is that both these state 
equations can exist at the same time; they refer to two different subsystems. 
The properties of the vacuum fluid are described entirely by the A term in 
the modified energy-momentum tensor T^ u in Eq. (2); as mentioned earlier 
it corresponds to a vacuum energy p vac = A/8nG and a vacuum pressure 
Pvac = —A/8nG. In addition to this, the matter/radiation fluid forms a 
separate subsystem. 

The point that we wish to stress is: although the present theory in itself 
does not make a definite prediction about the relative energy fraction residing 
in the matter component, some amount of matter has always to be present, 
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in order to support the magnetic field. In a pure de Sitter universe consisting 
of expanding vacuum (cosmological constant only) there would be no charge 
carriers from which a magnetic field could be generated. Even if we imagine 
that there was for some reason a magnetic field present at some instant of 
time, this field would disappear quickly due to the expansion. It gives no 
physical meaning to consider MHD in an empty, purely vacuum, universe. 

(5) It has recently come to our attention that Subramanian and Barrow 
have published a detailed study [H| of consequences of MHD in the early 



FRW Universe. They discuss various aspects of the damping problem and 
also the possibility of observational detection of the magnetic fluctuations. 
Cf. also the related Ref. [|(| . 
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A Remarks on the FRW case 

We will give a brief account of the analogous theory for a radiation dominated 
flat FRW universe, both to put our above results into perpective, and also 
because our results for the dispersion equations for the modes differ partly 
from those obtained in [|l(| . 



The line element still has the form (5), but the scale factor R now becomes 

v 1/4 

Po) t 1/2 - (Al) 



1811 

R(t) ^32vr y/\ 1/2 



Ro V 3 
The equation of state for the radiation dominated fluid is 

V = \p. (A2) 

The governing equations are the /i = equation (25), the /i = i equa- 
tion (26), and the Maxwell equation (34), as before. Taking Eq. (A2) into 
account, we obtain for the \i = equation 

i i a ( J R 4 p) + ^V-(pu) = E-J, (A3) 
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whereas the \i = i equation becomes 



4 -M R3 PUi) + \ d iP = (B k d k )B t - l -d t B 2 . (A4) 



3i? 3 uv ' l > 3 l ' v K > l 2 

Considering plane, linear, MHD waves on a uniform background in this uni- 
verse we now have, instead of Eq. (35), 

PoW«^, BSWoc^. (A5) 

We expand the wave quantities analogously to Eq. (41) (for simplicity 
keeping the same symbols as previously for the constants): 

«• = nje* , ^P=^ e , 5B* = ^ @ . (A6) 

We obtain the conservation equations 

^p R 4 k ■ u = Coo, (A7) 

4 C 

-(jjp R 3 u 0i - — -hi = (b • B )ki - (k • B )b 0i , (A8) 

o ort 

which together with Maxwell's equation (46) determine the dispersion rela- 
tions. We distinguish between the same cases as above: 

If Bo is perpendicular to k, then there is a longitudinal wave (uq2 = uqs = 
0) with dispersion relation 

^ft + S) k2 (A9) 

This equation, contrary to Eq. (47), describes a physical disturbance, since 
both phase velocity and group velocity are subluminal. 

If B is parallel to k, there is an ordinary longitudinal nondispersive 
luminal wave, as before. Moreover there is a transverse wave whose dispersion 
relation is 

cu 2 = -^-(k-B ) 2 . (A10) 

4po 
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Comparison between Eqs. (47) and (A9), and between Eqs. (52) and (A10), 
shows how the de Sitter- Zel'dovich case differs from the radiation dominated 
FRW case. Whereas the longitudinal wave (A9) agrees with Eq. (31) in [[IT]], 
there appears to be a deviation in case of the transverse wave, Eq. (A10). 
It is also noteworthy that Eq. (A9) is in agreement with the relationship 



viong = \Js 2 + v\, where now s = v/l/3 and va = y3/4p B is the FRW 
Alfven velocity corresponding to Eq. (A10). 
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